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Abstract 

An explicit convergent iterative solution for the lowest energy state of the Schroedinger 
equation with generalized A^-dimensional Sombrero-shaped potential is presented. The 
condition for the convergence of the iteration procedure and the dependence of the 
shape of the groundstate wave function on the parameters are discussed. 
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1. Introduction 



In this paper, the iterative method of Friedberg, Lee and Zhao[l] will be ap- 
plied to a generalized radially synnnnetric Sombrero-shaped potential in N- 
dimensional space with 

V(r) = \g\r'-rl)\r' + Arl), (1) 

where rg = (2 + A^)/3, and A are arbitrary constants. The same potential 
in the special case of = 1 has been studied recently[2]. As noted in Ref.[2] 
this investigation is stimulated by an interesting question raised by Roman 
Jackiw[3]. For arbitrary A^, the corresponding Schroedinger equation for the 
groundstate radial wave function is 

with k = {N — l)/2. The boundary conditions are 

V'(oo) = and ^''(O) = 0. (3) 

When g = 1 and A = 2 the solution of the groundstate has an analytical form 
as tjj{r) = 1^ with the eigenvalue = Tq. However, for arbitrary g and A 
the groundstate wave function has no analytical form. 

To apply the iterative method we introduce a trial function 0(r) satisfying 
another Schroedinger equation 

(-^|:^''|: + - Hrmir) = gEocf>{r) = {E- A)0(r), (4) 

where h{r) and A are the corrections of the potential and the groundstate 
energy. The explicit construction of the trial function 0(r) will be given in the 
next section. Here, we note some useful relations. 

Multiplying (4) on the left by ip{r) and (2) by 0(r), their difference gives 



Let 

^(r) = 0(r)/(r). (6) 
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The equations for /(r) and A can be deduced by using 

2k 1 2 df 



A. 

dr 



= 2r2^(/i - A)0V- (7) 



The integration of the left-hand side of (7) over r = to oo is zero. This 
gives the expression of the energy correction 

OO 

/ r'^^(lP'{r)h{r)f{r)dr 

A = ^ . (8) 

/ r'^^(j)'^{r)f{r)dr 

Introducing two iterative series {fni'f)} and {A^^} with n = 0, 1, • • •, we 
obtain the following iteration equations 

/ r^'^(l)^{r)h{r)fn~i{r)dr 

An = ^ , (9) 



Ur) = Un) - 2 / / x2V'(x)(A, - hix))fn-iix)dx, (10) 

where r^, could be chosen as r^, = or = oo. To ensure the convergency of 
the iterative method it is necessary to construct the trial function in such way 
that the perturbed potential h{r) is always positive (or negative) and finite 
everywhere. Specially, h{r) when r oo. In the following we construct 
the trial function for the iteration procedure. As we shall see, in the cases that 
we have examined the iterations give rapidly convergent results. 
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2. Trial Functions 



We first introduce a trial function 

0+(r) = f!l±^)'e-^'^°(^)-'^^W (11) 
satisfying the following Schroedinger equation 

and the boundary condition 

Moo) = 0. (13) 

Substituting (11) into (12), we connpare terms with the sanne power of g. From 
p^- terms we obtain 



<S^,(r) = ^/2^ = (r^ - rl)^r'^ ^ Arl (14) 

To ensure h+{r) satisfying the convergence condition, 5i(r) is defined in a 
special way to prevent terms with positive powers of r presenting in h+{r). 
For terms we have 

-1(2S'qS[ -S'^- 2ka^r'' + Ar^) - ka^ = Eq. (15) 

2 V / r=ro 

Introducing 

Eo = + 4'^ + 4'^ (16) 

and defining 

= -A;a2 (17) 

we write 

Since 5'o(ro) = we obtain 

4'^ = l^oiro) = rlVTTA, (18) 
= karoVTTA (19) 



and 



S[ = i^S'^ - Ei'^)/S', + [ka^fi^^^o - (20) 



Substituting <So(r) into (20) we have explicitly 



, . . _ r(3r^ + 2Arl - rg) - 2rl^/TTA^/r'^ + ^rg 
•^lin - 2(r2-rg)(r2 + ylr§) 



The expression for /i+(r) is 



(21) 



2^ ^ 2 (r + a)2 r(r + a) r(r + a) 



+fca,(r? - a^)^^^^ + ka'g (22) 
r(r + a) r(Vr2 + Ar^ + r) 

The condition 0'(O) = can be satisfied by introducing the trial function as 

(/){r) = 0+(r) + ^(j)-{r) for r < ro (23a) 

and 

0(r) = (l + ^0-(ro)/0+(ro))0+(r) for r > ro (236) 
where 0-(r) is defined as 

0_(r) = (^)'e-^'^°(-^)-^^W. (24) 

The parameter ^ is fixed to satisfy the condition 0'(O) = 0, namely 

0V(O)+^0-(O) = O. (25) 
Correspondingly 0(r) satisfies the Schroedinger equation (4) with 

h{r) = h+{r) for r > vq (26) 

and 

2 2 

/i(r) = h+{r) + 2p^fEo + ka " ~ + Arg 

V r(r + aj 



-ka'—r^=^ -U-ir)/4>ir) (27) 



for r < tq. It is interesting to notice that the condition (25) for 0'(O) = also 
ensures h{r) to be finite when r ^ 0, which is necessary for the convergency 
of the iteration procedure. 

By integrating (14) and (21) we obtain So{r) and *Si(r) as 



^o{r) = o^yr^ + Ar^{2r^ + Ar^ - 4rl) - -{A^Vq + Mrg) ln(r +\/r^ + Ar^) 
8 8 

(28) 

Q f \ w I \ ^ ^ 2 , A 2^ , A ka Vi-^A^r^ + Ar^ + r + Atq 
5,(r) =l„(r+r„)+-ln(r +Ar,)+i-+-) In + _ , + 

(29) 

Substituting them into (12), (24) and (23a-b) gives the final expression of the 
trial function. From (21) we can also reach 

^(<5f-<sr) = --^-i,--^+ 



8(r2 + Ar§)2(a + /3) 2r(r + a){r^ + Ar^) a' + /3' 



2(^2 ^ Ar§)(^r2 + ^r§ + roy/TTAy 
kar{2yjr'^ + Arg + ro\/l + ^) 



(r2 + ylr§)3/2(^r2 + ^rg + roVTTAy 



(30) 



where 



7 = — V = - ^ (31) 

7 (r2-r§)2' ^ (r2-r§)2 



with 

a = 15r^ + (18A - 6)rV^ + (8^^ + 12A + 7)rVJ + (8^^ + 2A)r[ 

(3 = 8Vl + ^r^r(3r2 + (2^ - l)r^yr2 + ^r§, 
a' = r(3r^ + (2A - l)r^) 

and 
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Explicitly 

7 = 225r^ + 270(1 + 2A)rV^ + 3(188^2 + 216A - 5)rV^ 
+36^(8^2 + 10^ - l)rV§ + 4(16^^ + 8A^ + l)rg 

and 

y = 9r^ + 3(4A - l)rV^ + 4^(1 + A)^. 

Substituting (30) into (22) and (27) gives the final expression of h. With 
above results for h and we are ready to perform the iteration. 

3. Numerical Result 

Starting from the above defined trial function 0(r) and the related h{r), we 
can perform the iteration based on (9) and (10). Our numerical results show 
that the finally obtained wave functions and eigenvalues for the groundstate 
convergent nicely. Let us take N = 3 bs an example. 

For g = 1 and A = 2, the exact solution of the groundstate is 

with Eq = Tq. Starting from the trial function defined by (23a) and (23b), 
with arbitrarily chosen parameter a and the corresponding parameter ^ fixed by 
(25), the iteration procedure (9) and (10) gives the final convergent result of 
the wave function and the eigenvalue of the groundstate, which is consistent 
to the exact solution. The trial function and the convergent groundstate wave 
function after the iteration are plotted in Fig. 1. It is interesting to observe the 
transition of the shape of the wave function for the trial function with maxima 
at a finite r to the final convergent one with only one maximum at r = 0, as 
the exact groundstate wave function should be. This answered the question 
raised by R. Jackiw[3] in A^-dimensional case: Even the trial function proposed 
has its maxima at r > the iteration procedure would still reach the exact 
solution of the groundstate wave function with its only maximum at r = 0. 

As examples the obtained groundstate wave functions after the iteration 
procedure are plotted in Figs. 2 and 3 for ^ = 2 and g = 0.5, 1 or 2, and for 
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g = 1 and A = 1, 2 or 3, respectively. In Table 1 the corresponding eigenvalues 
of the groundstate obtained fronn the iteration are listed for different paranneters 
g and A. It can be seen that the iterative energy series convergent quite nicely. 
From the figures it is interesting to see the transition of the form of the obtained 
groundstate wave functions from the shape with maximum at r = to the 
one with maxima at a finite r, becoming a degenerate groundstate, when g 
increases from < 1, passing 1 to > 1 for A = 2, or when A increases from < 2, 
passing 2 to > 2 for p = 1. The results seem to show that the groundstate 
wave functions in the region g < 1 and A < 2 have the shape with only one 
maximum at r = 0, while in the region outside the wave functions become 
degenerate at a finite r. Their maxima move to larger r when the parameters 
g and A increase further. 
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Table 1. Eigenvalues of groundstates in the iteration for = 3 



g A 


gEo 


gEi 


gE2 


gEs 


gE, 


gE, 


0.5 2 


-0.4300 


1.3963 


1.3795 


1.3775 


1.3773 


1.3772 


1 2 


-8.6479 


2.1523 


2.1516 


2.1517 


2.1517 


2.1517 


2 2 


5.5581 


4.1362 


4.0976 


4.1108 


4.1092 


4.1094 


1 1 


-2.3537 


1.8920 


1.8473 


1.8402 


1.8393 


1.8392 


1 3 


3.6773 


2.4675 


2.4353 


2.4425 


2.4417 


2.4418 
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Fig. 1 Trial Function 0(r) and Groundstate Wave Function ^(r) 
for = 3, ^ = 1 and A = 2. 




Fig. 2 Groundstate Wave Function ^(r) for = 3, gi = 1 
and A = 1 (thin), 2 (middle) and 3 (thick). 
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Fig. 3 Groundstate Wave Function ^(r) for N — 3, A — 2 
and g = 0.5 (thin), 1 (middle) and 2 (thick). 
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